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Abstract 

We propose a scheme for the quantum non-demolition (QND) measurement of Fock states of a nanomechanical resonator 
via feedback control of a coupled circuit QED system. A Cooper pair box (CPB) is coupled to both the nanoresonator and 
microwave cavity. The CPB is read-out via homodyne detection on the cavity and feedback control is used to effect a non- 
dissipative measurement of the CPB. This realizes an indirect QND measurement of the nanoresonator via a second-order 
coupling of the CPB to the nanoresonator number operator. The phonon number of the Fock state may be determined by 
integrating the stochastic master equation derived, or by processing of the measurement signal. 
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I. INTRODUCTION 

An important benchmark for quantum control is the 
ability to prepare and detect a harmonic oscillator in 
a Fock state, an energy cigenstatc. This may soon be 
achieved in nanomechanical resonators^, which are more 
usually prepared in quasi-classical displaced thermal 
states due to the ubiquity of environmental interactions 2 -. 
The ability to prepare a mechanical resonator in a Fock 
state will provide access to processes that mark the 
quantum-to-classical transition^. 

Fock states have already been prepared in other 
instances of harmonic oscillators. Using a cavity 
QED system, few-photon Fock states of an electro- 
magnetic mode have been prepared deterministically 4 , 
via state reduction^ and via quantum non-demolition 
measurement^. Few-photon optical Fock states have 
also been generated conditionally using parametric down 
conversion^. Taking advantage of the higher level of 
control afforded by circuit QED, Fock states contain- 
ing up to fifteen photons have also been deterministi- 
cally generated ^ 10 ! 11 . The motional state of a trapped 
ion provides another instance of a harmonic oscillator, 
and indeed Fock states have been prepared by tuning 
the interaction between the ion's motional state and its 
internal electronic levels via an external laser—. 

One may consider analogous schemes for the genera- 
tion of such states of macroscopic mechanical resonators. 
Indeed, in an extraordinary experiment, Fock states of 
a mechanical oscillator, in the form of a suspended film 
bulk acoustic resonator resonantly coupled to a supercon- 
ducting phase qubit, have been demonstrated^. How- 
ever, it would still be desirable to generate mechanical 
Fock states in a lower-frequency mechanical resonator 
with a clear separation of mechanical and electromag- 
netic degrees of freedom. In order to conditionally pre- 
pare a Fock state of such a resonator, one requires a cou- 
pling to its number operator. In general, this is difficult 
to realize, though an optomechanical "membrane-in-the- 
middle" system does so^. In other key ingredient 

is the coupling of the resonator to a qubit. Coupling be- 
tween a nanoresonator and a superconducting qubit has 



recently been demonstrated by LaHaye et a/.—. Irish and 
Schwab proposed a scheme for the generation and detec- 
tion of a Fock state of a nanoresonator via coupling to 
a Cooper pair bo»i£. Jacobs et al. have considered the 
continuous measurement of a nanoresonator using cir- 
cuit QEDii; the nanoresonator is coupled to the qubit, 
and the qubit is read out using a superconducting mi- 
crowave cavity^. Here the linear coupling betweeen the 
nanoresonator and the qubit is brought into resonance 
by driving the nanoresonator at a frequency equal to the 
detuning between the nanoresonator and the qubit. It 
is shown that, in spite of the linear coupling, the effect 
of the measurement is to localize the nanoresonator onto 
its Fock states. A dispersive coupling to a qubit with the 
nanoresonator number operator has also been considered, 
with the proposed read-out in this case being through a 
single-electron transistor or a quantum point contact^. 
The measurement of Fock states by direct coupling to 
a single-electron transistor has also been investigated^. 
Alternative proposals for conditional Fock state gener- 
ation in nanoresonators are based on engineered cross- 
Kerr interaction ;) 21 ' 22 , in analogy with the situation com- 
monly encountered in an optical setting 2 ^. 

The system under consideration here is composed of a 
nanoresonator, a superconducting microwave cavity and 
a Cooper pair box. Both the nanoresonator and the mi- 
crowave cavity are separately coupled to the Cooper pair 
box (CPB)2i. The direct coupling between the nanores- 
onator and the cavity could also be considered, and it 
may be shown that resolved sideband is achievable us- 
ing a similar analysis 2 ^. There is a resonant second- 
order coupling between the nanoresonator and the CPB, 
which gives a Hamiltonian coupling between the nanores- 
onator number operator and the qubit inversion o~ z (in 
the energy eigenbasis of the qubit) as Wba z . In the adia- 
batic limit, the qubit coherence (<r_) may be continuously 
read-out via homodyne detection on the output of the 
microwave cavity, and via feedback control of the CPB 
charge bias a non-dissipative measurement of o~ v may be 
created 2 ^. Such feedback has been considered before in 
the context of stabilizing pure states of a qubilj 2 ^. Thus 
we have an indirect measurement of the nanoresonator 
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phonon number, and may conditionally prepare a Fock 
state of the nanoresonator—. In principle, the condi- 
tional evolution of the system may be reconstructed by 
integrating the stochastic master equation, though the 
phonon number may be extracted (and quantum jumps 
observed) from the oscillations of a y by filtering the ho- 
modyne measurement record. Note that feedback con- 
trol of nanomcchanical systems has previously been con- 
sidered for the purpose of cooling the mechanical mode 
toward its ground stat e 28 ' 29 and for squeezing of its quan- 
tum fluctuation o 30 ' 31 . 



II. SYSTEM AND HAMILTONIAN 

The Hamiltonian describing this system, in the 
Schrodinger picture and in the charge basis of the CPB, 
is 

H s = frw c a)a + hwmtfb + y a z + 

+hg (a + a^) a z + h\(b + tf)a z (1) 

where uj c is the resonance frequency of the cavity, uj m is 
the resonance frequency of the nanoresonator, and the 
CPB is described by he = -4£ c (l - 2n g ) and HA = 
-£ , ,/cos(7T$ e /$ )- Here E c = e 2 /2Cs is the charging 
energy of the box, n g = n™ + n g} n™ = C g n V g n /2e is 
the dimcnsionless gate charge due to the nanoresonator 
voltage, n c g = C g V g /2e is the dimcnsionless gate charge 
due to an additional control voltage, Ej = I c $>q/2tt is 
the Josephson energy of each of two Josephson junc- 
tions between the CPB and its superconducting reser- 
voir, 4> e is the external flux threaded through the two 
Josephson junctions, <£>o is the magnetic flux quantum 
and I c is the critical current of each Josephson junction 
to the CPB. The total capacitance of the box is given by 
Cs = 2C,j + C g + C g n + C g where Cj is the capacitance of 
each Josephson junction, C g is the coupling between the 
box and the central conductor of the microwave cavity, 
is the capacitance between the nanoresonator and 
the box, and C g is the capacitance between the control 
gate and the box. 

It is assumed that the CPB is configured as a charge 
qubit, meaning that Ec >> Ej. The couplings between 
the CPB and the cavity 3 ^, and between the CPB and the 
nanoresonator— are given by, respectively, 

s '«= e ty¥< P) 

Here c is the capacitance per unit length of the microwave 
cavity, L is the length of the microwave cavity, Aa; is the 
half-width of the nanoresonator ground-state wavefunc- 
tion, and d is the equilibrium separation of the nanores- 
onator and the CPB. 

First we rotate into the energy eigenbasis of the un- 
coupled qubit; the natural basis for a treatment of dis- 
sipation. Additionally transforming into an interaction 



picture with respect to Hi = huj c a) a, 

+ Ml , x N / e A \ 

H h = huj m tfb + —a z + h\ (b + W) f —a, - —a x J 

+hg (ae-^ + a)e^) (j^a z - ^ , (3) 

where the CPB energy splitting is ft = \J t 2 + A 2 . 

In the absence of resonant terms in the coupling 
between the nanoresonator and the CPB, we make a 
Schrieffcr- Wolff transformation 3 ^ to extract a resonant 
second-order coupling; Hf™ = S j( Hj 1 S where S = 
exp [i(\A/n 2 )a y (b + . This is essentially equivalent 
to considering the dispersive limit of the coupling be- 
tween a qubit and an oscillator. Retaining terms to lin- 
ear order in AA/J7 2 , transforming to another interaction 
picture with respect to H2 = h(u> c /2) a z + ku! m Wb, and 
employing the rotating-wave approximation to neglect 
rapidly oscillating terms, we find 

H l2 = h(5 + xtfb)<T z + hg' (atr + + a t (T_) , (4) 

where 5 = (£1 — cj c ) /2 is the detuning between the CPB 
and the cavity, g' = —gA/Q is the effective coupling to 
the cavity, and x = 4A 2 A 2 /f2 3 is the effective coupling to 
the nanoresonator. The cavity has a Jaynes-Cummings 
coupling to the CPB, which itself is dispcrsively coupled 
to the nanoresonator. 



III. CONDITIONAL DYNAMICS UNDER HO- 
MODYNE DETECTION 

We consider the conditional evolution of the system de- 
scribed by Eq. (01 under homodyne detection of the field 
output from the microwave cavity. There exist proposals 
for near quantum-limited measurement of the quadra- 
tures of a microwave fiel d 35 ' 36 . It is assumed that the 
microwave field, leaking out of the cavity at a rate [i, is 
subject to homodyne detection with respect to a local os- 
cillator phase 9 and with an efficiency 77, that the nanores- 
onator is damped at a rate 7 into a thermal environment 
characterised by the occupation n^, and that the CPB 
experiences spontaneous emission at a rate T q . The ef- 
fect of the Schrieffcr- Wolff transformation on dissipativc 
terms may be neglected under the joint assumptions of 
weak damping and AA/f2 2 << I. The master equation 
describing the conditional evolution of the system is 

dp = ~[H l2 ,p]dt + 'y(n°+l)V[b]pdt 

+jn° m V [6 t ] pdt + ptD [a] pdt + T q V [tr_] pdt 
+ ^qpU[ae- l9 ]p dW, (5) 

with dr = (ae~ l6 + a^e %e ) dt + dW/y/rfp, as the measure- 
ment record increment (where dW is the Wiener incre- 
ment), T>(s)p = sps^ — S p _ ips^ s as the dissipativc 
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superoperator and H[s]p = sp + ps' — tr (sp + ps^ p as 
the measurement superoperator—. 

In the good-measurement limit, the cavity damping 
rate is large, p » 5,x,g', and we may adiabatically 
eliminate the cavity mode. This is easily done by solv- 
ing the quantum Langevin equation in the limit of large 
damping to find a = (—2ig'/p) er_ . Substitution into the 
conditional master equation above leads to a conditional 
master equation for the reduced density operator describ- 
ing the CPB and the nanoresonator, p r (t) = tr cav i ty p(t), 



-i [(5 + x^b)a z ,p r ] dt + r q T) [a-] p r dt 
+7 « + l)V [b] Pr dt + in ° m V [fct] Pr dt 

+TV [cr_] p r dt + ^/rjrn [cr_ e -* (e+7r/2) ] p r d\f6) 



dp r = 



where the measurement rate is T = 4g' 2 / p, and the mea- 
surement record increment (setting 6 = — n) is 



dr 



(a y ) dt 



dW 



'//' 



(7) 



The associated measurement current is defined through 
dr = I(t) dt. 



then according to the quantum theory of continuous 
feedback^, the system evolution is given by 



da = ~[H,a}dt-^\[jF + Fc,a]dt 

A 2 1 - ?/ 
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V [F] adt 

(11) 



Now applying the Hamiltonian feedback of Eq. (fT0| to 
the system described by Eq. © with the identification 
c = i<7_ in Eq. pip , we find 

dp r = —i [((5 + \b b) &z> Pr] dt + T q T> [aJ\ p r dt 
+7 (n° n + l) V [b] p r dt + jn° n V [tf] p r dt 
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x - 
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V [F] p r dt + V 



Vr . 



-H 



iy/r)Ta- - iXF/y/ip? p r dW. 



p r dt 



(12) 



In general, the feedback will alter the Hamiltonian, dis- 
sipativc and measurement dynamics. 



IV. FEEDBACK CONTROL 



V. NON-DISSIPATIVE MEASUREMENT 



Now Eq. © describes a dissipative measurement of 
the CPB; the measurement perturbs the observable from 
which we wish to extract information about the nanores- 
onator. One may consider using state-based feedback 
control^ of the CPB to optimally extract information 
from the measurement; for example, one may lock the 
net detuning (S + xb^b) to zero such that the CPB ob- 
servable with a reactive response to the detuning exhibits 
optimal sensitivity to the nanoresonator state. 

However, a simpler alternative is afforded by Hamil- 
tonian feedback control using the measurement record^. 
Assuming that the feedback may be applied fast on the 
time-scale of system dynamics, we may apply the Marko- 
vian homodyne feedback master equation. Consider a 
density matrix a evolving according to the homodyne 
measurement master equation in Lindblad form, 



da = - - [H, a] dt + kV [c] adt + ^i]kU [c] a dW, (8) 

and having the associated homodyne measurement cur- 
rent, 



I(t) = (c + J)(t)+at)/Vvk, 



(9) 



where k is the measurement rate and £(i) = dW/dt. 
Then, if one applies Markovian feedback using the ho- 
modyne measurement current according to the Hamilto- 
nian, 



Now it is possible, in the case of unit detection 
efficiency, to create a non-dissipative measurement of 
<j y from a dissipative measurement 2 ^. In order to 
achieve such a measurement here, we require that XF = 
(r)T/2)a x . The corresponding feedback Hamiltonian, 
given by Eq. (|10j) , is in a rotated qubit basis and an inter- 
action picture, as per Eq. In the Schrodinger picture 
and charge basis of Eq. ([1]), the required feedback term 
is -h£ fb (t)cosClta z where £ fb {t) = (r]TQ/A) I(t). That 
is, the homodyne measurement record must be used to 
modulate the CPB charge bias, controlled through Vg, 
to create a non-dissipative measurement of a y . 

The associated homodyne feedback master equation is 

dp r = —i [(5 + xb^b) a z ,p r ] dt 
+7 (n° m + 1)V [b] Pr dt 



T q V [a 



p r dt 
jn° m V [6+] p r dt 



-TV 



V 
2 



P7? 

rdt + (1 — 17) [Cx] Prdt 



+ ^/t 1 TH [a y /2} p r dW. 



(13) 



H fb (t) = HXFI(t), 



(10) 



Here there is no alteration to the Hamiltonian dynamics, 
though for inefficient detection, the feedback also intro- 
duces a dephasing of a x , and the measurement becomes 
dissipative. 



VI. SIMULATIONS 

It is assumed that the nanoresonator has a resonance 
frequency o; m /27r = 10MHz and a mass m = 10~ 15 kg, 
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FIG. 1: Evolution of: (a) phonon number (fo^fe) and (b) 

phonon number variance Var (b^fe) = ((Vfo) 2 ) — {b^b} 2 for a 
typical trajectory assuming unit detection efficiency (rj — 1). 
Parameters are as given in the body of the text. The inset 
shows the conditioning onto a Fock state that occurs on a 
time-scale short compared to that required for the observa- 
tion of retroactive quantum jumps. 



corresponding to ground-state half-width of Ax = 29fm . 
These are reasonable parameters for a nanoresonator 
fabricated within a superconducting microwave cavity^ 4 -. 



The CPB is characterized by A = 
such that e = 0, and hence f2 



5 x 10 10 s 



1/2 



5 x 10 iu s 



nant with the cavity at ui c = 5 x 10 s ; such values 
are typical for circuit QED experiments^ 8 -. The cavity 
damping rate is fj, = 10 s , and the assumed mechan- 
ical quality factor is, an admittedly high, 2 x 10 7 . The 
bare qubit-nanoresonator coupling is X/2tt = 1.5MHz, 
and the bare qubit-cavity coupling is <?/27r = 200kHz, 
both of which are achievabl e 15 ' 18 . The effective CPB- 
nanoresonator coupling is \ = 2-56 x 10 3 s _1 , and the 
effective qubit-cavity coupling is g' = —7.56 x 10 5 s _1 . 
The effect of the cavity is quantified by the measurement 
rate T = 2.29 x 10 5 s _1 , and spontaneous emission from 
the qubit occurs at a rate r g = 10 4 s _1 . 

Simulations of the mixed-state evolution of the 
nanoresonator-qubit system according to Eq. (|13[) were 
performed using an explicit Milstcin scheme 41 -. Since 
this is an instance of a partially-observed system, con- 
ventional pure-state unravellings 4 ^ of the master equa- 
tion are not applicable, though techniques for pure-state 
unravellings using the concept of ostensible trajectories 
have been developed 4 ^. 

Fig. 1(a) and Fig. 1(b) show, respectively, the evolution 
of (b^b) and Var (b^b) over the course of a typical quan- 
tum trajectory assuming unit detection efficiency. The 
initial state in this case is a thermal state with fi = 2 and 
no coherent amplitude. This is close to the level of cool- 
ing recently achieved in a microwave cavity optomechan- 
ical system through resolved sideband cooling 4 ^. From 
Fig. 1(a) it is seen that the effect of the measurement 
is to condition the nanoresonator onto a state with an 
integer (tfb), and from Fig. 1(b) it is seen that Var (&^&) 
goes to zero, verifying that the nanoresonator has been 
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FIG. 2: The interaction picture evolution of the qubit observ- 
able (a y ) corresponding to a short interval of the simulation 
shown in Fig. 1. The dynamics are the result of a compromise 
between projection onto the eigenstates of a y and oscillation 
due to the free Hamitonian, determined by the relative mag- 
nitude of T and xi^b). 



projected onto a Fock state by virtue of the indirect mea- 
surement. In the absence of the measurement, the vari- 
ance would tend to (n^ + l), in equilibrium with the 
thermal mechanical heat bath. After some time, jumps 
between adjacent Fock states are observed, correspond- 
ing to the emission and absorption of quanta to/from 
the thermal mechanical heat bath. These have previ- 
ously been termed retroactive quantum jumps 4 ^, and 
arise from the interplay of system and measurement dy- 
namics, in spite of the lack of explicit jump terms in the 
stochastic master equation. 

Fig. 2 shows the interaction picture evolution of the 
measured observable a y over a short interval of the simu- 
lation shown in Fig. 1. The measurement acts to put the 
qubit into an eigenstate of a y , while the free Hamiltonian 
a z component drives oscillations of (a y ) at a frequency 
dependent on the phonon number of the nanoresonator. 
The dominant effect is determined by the relative mag- 
nitudes of the measurement rate T and the Hamiltonian 
coupling coefficient x(tfb}. In Fig. 2, both effects are 
observable, and it is clear that by monitoring (cr y ), one 
acquires information about (b^b). 
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FIG. 3: Comparison of phonon number variance for typical 
trajectories following projection onto a Fock state, for efficient 
and inefficient detection. Greater fluctuations are observed in 
the case of inefficient detection. 

In the presence of detector inefficiency, one still ob- 
serves projection onto the nanoresonator's Fock states, 
though a longer time is required for this to be achieved. 
As shown in Fig. 3, the fluctuations in phonon number 
variance observed in the case of inefficient detection tend 
to be greater than those observed in the case of efficient 
detection. 
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VII. CONCLUSIONS 

Microwave cavities provide a natural environment for 
the application of quantum control protocols; the emis- 
sion of a qubit can be effectively confined by the mi- 
crowave cavity, and one can feed back onto the qubit 
or the microwave cavity on a time-scale comparable to 
the system frequencies. Here we propose a system that 



exploits this possibility by using feedback control of a 
Cooper pair box qubit to implement a quantum non- 
demolition of the Fock state of a nanomcchanical res- 
onator. The most significant experimental challenges are 
the fabrication of mechanical resonators with sufficiently 
high quality factors, and the efficient measurement of a 
quadrature of the microwave cavity field. 
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